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In this work we classify the singularities obtained from the Gibbs potential of 
a lattice gas model with three components, two order parameters and five control 
parameters applying the general theorems provided by Catastrophe Theory. In par- 
ticular, we clearly establish the existence of Landau potentials in two variables or, in 
other words, corank 2 canonical forms that are associated to the hyperbolic umbilic, 
-D+4, its dual the elliptic umbilic, -D-4, and the parabolic umbilic, D5, catastrophes. 
The transversality of the potential with two order parameters is explicitely shown for 
each case. Thus we complete the Catastrophe Theory analysis of the three-component 
lattice model, initiated in a previous paper jjj]. 
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I. INTRODUCTION 

The importance of phase transitions with several order parameters is very well known in different 
branches of physics Q . A great amount of theoretical work has been done in order to understand and 
construct accurate phase diagrams. As is well known, two different approaches are usually employed, 
one more phenomenological by introducing Landau polynomial potentials which try to describe ex- 
perimental singular behaviors and the second one applying Catastrophe Theory (CT) (also known 
as Singularity Theory and thus adopting a more methodological point of view. Even in the 

second case there are several ways to deal with phase transitions and diagrams. Most works adopting 
the second point of view start with the canonical unfoldings as given and base their treatments on 
the effect of perturbations leading to preserving or not the internal symmetry of the system consid- 
ered or, in other words, they focus their analysis in the symmetry-breaking character of some phase 
transitions. The procedure we adopt here is different since we begin with a thermodynamic potential 
(for example, the Gibbs potential), assuming a mean field approach, and we apply an algorithm or 
program according to the general mathematical theorems established by CT, in order to extract all 
the topological information of the original thermodynamics potential. In a recent paper we have ap- 
plied this CT program to a three component model phase diagram and we have found for a one order 
parameter potential the highest singularity with codimcnsion five to be the wigwam or Aq catastrophe 

0- 

Our CT program could be briefly stated as follows. Let H(x,ir) be a real function with state 
variables x\, . . . ,x n (x £ R") and control parameters 7Ti, . . . , n r (n £ K r ); that is, H : H" +r — » K. 
Then we must take the subsequent steps: 

1. We pick (xo,7To) such that £0 is a degenerate critical point of H(x,ttq) and we consider the 
unfolding h(x,n) = H(x + xq,tt + ttq) — H(xo,ttq) and g(x) — h(x,0) to translate x$ to the 
origin of coordinates. 

2. We calculate the determinacy and codimcnsion of g from the fc-jet of g. Of course, if g is k- 
determinate then g ~ j k (g), that is, the function g is equal to j k (g) up to a change of coordinates. 
In other words, they are equivalent and have qualitatively the same properties; therefore, both 
have the same codimension, cod(g) — cod(j k (g)). 
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3. We study the /c-transversality of h. If this function is fc-tranversal we can affirm that h and 
the canonical form of the unfolding of g are isomorphic. Then we can replace the original H 
function by this canonical unfolding. If not, we can state that the H function is not susceptible 
to be studied by CT. 



CT has usually not been applied in a rigorous way using all these concepts and theorems needed 
for its correct implementation. This Catastrophe Program proposed here provides a very useful and 
systematic way to examine and classify with not very much effort general behaviors of physical systems. 
In particular, we emphasize the study of transversality of the actual thermodynamical potentials which 
guarantees that those simple forms (polynomial potentials or canonical forms) represent indeed up to 
a diffeomorphism the original thermodynamical potential. Following this CT program we do not need 
to invoke any convention (for example, delay or Maxwell convention) in order to classify degenerate 
or non-degenerate critical points on the state variables space. Both conventions are not intrinsic to 
CT. Only when we deal with the time evolution or when dynamical considerations about the physical 
system are considered, a given convention could be advisable. In particular, when the order of a phase 
transition needs to be determined a convention is necessary because the transition occurs when an 
appropriate separatrix in the control parameter space is crossed. 

Here we will focus our attention on the lattice gas model for a system with three components 
which simulates, in particular, a binary fluid mixture. A wide literature has already been devoted to 
it from different points of view [fz| — p~0|, and references therein], restricted to the case with one order 
parameter. In fact, very few studies with these methods can be found for phase transitions with two 
order parameters plj . This is rather surprising since the Landau potential for the three-state Potts 
model, which is a particularly important three-component model, has long been known to have two 
order parameters . We shall perform a complete CT study of the case with two order parameters, 
which is the maximum number for this model. 

This work is organized as follows: In section 2 we describe the thermodynamical potential to be 
analyzed, give its physical interpretation and discuss general stability questions which help connect 
usual concepts in Thermodynamics with those in Catastrophe Theory. In section 3 we apply the CT 
program to the potential previously introduced by considering the singularities with corank equal to 
2 and establish the elementary catastrophes associated. In Section 4, we analyze in more detail the 
Potts model as a particular case. The last section is devoted to a discussion of the previous results 
and of the structure of the phase diagram entailed by them. 



II. DESCRIPTION OF THE GIBBS POTENTIAL 



In the mean field theory, the Gibbs potential is a function of the concentration of two of the 
three components and depends on three thermodynamical parameters, which can be taken as the 
temperature and the chemical potentials of the two components, and on three molecular parameters. 
The phase diagram deduced from this function is an accurate description of the system, except close 
to the (multi)critical points, where fluctuations become important and alter significantly the mean 
field theory predictions. For this reason, the Gibbs potential has been the basis for determining the 
overall phase diagram |^|| . 

According to Ref. J8| a phenomenological model for a ternary mixture is obtained by assuming the 
Gibbs potential in the form 

G = N[a'yz + b'xz + c'xy + RT {x\nx + ylny + z\nz)], (1) 

where N — N x + N y + N z gives the number of total moles and N x , N y and N z the moles of each 
component. The variables x,y,z are the mole fractions defined by x = N x /N, y — N y /N and 
z = N z /N; and hence we have the constraint 

x + y + z = 1, with < x, y, z < 1. (2) 
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This model can be derived from the mean field theory of a lattice model Hamiltonian with variables 
taking three different states, representing the molecules of the three components (T^|. Then a',b' 
and c' represent molecular interaction parameters. The two order parameters to be considered in the 
sequel can be roughly thought of as the two independent differences between concentrations allowed 
by the constraint (|J). 

Several systems of interest are described by this Gibbs potential, Eq. (|l|): A ternary mixture at 

constant volumen, for example, a mixture of metals; a spin lattice where the molecules have spin one; 

a binary fluid mixture, where one of the three states represents a vacancy instead of a new molecule 

and the corresponding concentration is associated to a variable total volumen. In the last case, the 

possible phases are vapor, misciblc liquid mixture and inmiscible liquid mixture. The convenient 

extensive variables are the specific volume v and the relative concentration x — — r— of the two fluids 

. . x+v 

and the intensive variables are the pressure and the chemical potential of one of the fluids. Moreover, 

the thermodynamical potential Eq. (|l|) depends on T, v and x and is therefore the Hclmholtz potential 

F(T,v,x). This system is perhaps the most interesting for applications, given the great amount of 

experimental data on binary fluid mixtures [ pT|Jl5| . However, the potential (jl]) is not the most popular 

for fitting data; a related form which has similar dependence on the relative concentration of the two 

fluids but is of Van der Waals type for the volume is usually considered instead. We believe that this 

form, which is much more difficult to analyze, gives essentially the same qualitative behavior. 



III. ANALYSIS OF THE GIBBS POTENTIAL WITHIN THE FRAMEWORK OF 

CATASTROPHE THEORY 



Let us consider the reduced form of the Gibbs potential obtained from Eq. (|l|) and dividing by 
NRT, 

G(x, y, z, a, b, c) = a yz + bxz + cxy + xhix + yhiy + z In z (3) 

where now the new parameters a,b,c are defined with respect to the old ones a', 6', c' dividing them 
by RT. From the constraint Eq.(^|), we build a new function H of two variables such that 

H(x, y, a,b, c) = a y{\ — x — y) + bx(l — x — y) + cxy + xhix + yhiy 

+ (l-x-y)ln(l-x-y). (4) 



The mean field theory prescription is then to minimize the non- equilibrium Gibbs potential H — 
HxX~ y, 

H(x, y, a, b, c, [i x , [i y ) = ay(l — x — y) + bx(l — x — y) + cxy + xlnx + y lay 

+ (1 - x - y) ln(l - x - y) - fj, x x - y, y y, (5) 

with respect to x and y, where [i x and [i y are related to differences between the chemical potentials 
of the three components ||. 

CT will be applied to the H(x,y,a,b,c, /i x , /J, y ) function to classify the corank-2 singularities, at 
the generic point (xo,?/o) which moves on the triangle xq > 0, yo > and 1 — xq — yo > 0. CT 
conventionally uses the origin of coordinates as the point where singularities occur. Therefore, we 
shall translate the function H in order to have the singularities at the origin. This translated function 
is written now as 

h(x,y,a,b,c,fx x ,fx y ) = H(x + x ,y + y ,a + a ,b + b Q ,c + c , fi x + fi X fi,^ y + /x^o) 

- H(x ,yo,ao,b ,co,fJ, x ,o,fiy,o) (6) 

and the germ of the unfolding h is 

g(x,y)=h(x,y, 0,0, 0,0,0). (7) 



3 



Now according to the CT program proposed in our previous paper jjj and mentioned in the Introduc- 
tion, the following steps will be developed for our problem: 

1. Find conditions for which the origin of coordinates is a degenerate critical point of corank equal 
to 2 of g. 

2. Show that g is equivalent to p (polynomial), shortly expressed as g ~ p, which means that there 
is a diffcomorphism (p such that g = p ■ tp. The Dcterminacy and codimcnsion of g will also be 
evaluated. 

3. Find a canonical unfolding p (polynomial) of p by means of a basis of the vector space (x, y) / A(p) , 
which can be understood as the space of perturbations (see appendix B). 

4. Find a canonical 5-unfolding p\ of g. 

5. Establish the fc-transversality of h and if it holds, then by using the main theorems about 
fc-transversality, h and p\ will be isomorphic. 



A. Condition for a degenerate critical point at the origin. 

Now the first step is to write the conditions for which the origin of coordinates (0,0) is a degenerate 
critical point with corank equal to 2 of g. This can be done by equating all the first and second partial 
derivatives of g to zero at the point (0,0) (hypothesis So) — thus the Hessian of the germ g vanishes as 
well. This leads to five conditions among the variables and parameters xq, yo, ao> fro, Co, ^x,o and fi y ,o- 
From that system of five equations the following relations can be extracted 

«o = 2 ^o 1 +Vo *)> 
b o = 2 ( z o 1+x o 1 )^ 

co = \(x 1 +y 1 ), (8) 
x Zg 1 + zqx^ 1 +y Q Xq 1 ) + liix -lnz , 

XqZq 1 +z y ' 1 +x y^ 1 ) +lny -lnz , (9) 

with the definition zq — 1 — xq — yo > 0. Eqs. (||) give rise to a surface with parameters xq and yo 
fulfilling the condition for a degenerate critical point with corank 2 at (0,0) (Fig. 2). 



B. Classification of the germ g. 

With the hypothesis So, the 3-jet of g (Taylor expansion truncated beyond terms of degree 3) 
around the point (0, 0) is 

f{g) = ^x 3 (-Xo 2 + z 2 ) + ^x 2 yz^ + ^xy 2 z- 2 + ^y 3 (-y- 2 + z 2 ) 7 (10) 

This 3-jet can be considered now as a homogenous polynomial of degree equal to 3 and can be rewritten 
as 

j 3 (g) = aix 3 + a 2 x 2 y + a 3 xy 2 + a 4 y 3 , (11) 



and 



Mz,o = ^(-yoz Q 1 
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with ai = (1/6) (— x^ 2 + Zq 2 ), a 2 — a 3 = (1/2) z^ 2 and a 4 = (1/6) {—y^ 2 + Zq 2 ). We know from 
Singularity Theory that a general non-null homogeneous polynomial of degree 3 is equivalent by a 
linear transformation to one and only one of the following germs: x 3 — xy 2 , x 3 + xy 2 , x 2 y and x 3 
]l6| , p7t . We must remark that these x and y are not to be identified with the initial physical variables 
though they are linearly related to them. The application of this lemma to Eq. (|ll| ) to classify the 
3-jet is given in Appendix A and here only the final conclusions will be summarized: 

1. If (1 - 2x ) (1 - 2y ) (1 - 2z ) = then f{g) ~ x 2 y. 

2. If (1 - 2x ) (1 - 2y ) (1 - 2z ) > then j 3 (g) ~ x 3 + xy 2 . 

3. If (1 - 2x ) (1 - 2y ) (1 - 2z ) < then j 3 (g) ~ x 3 - xy 2 . 

All of these three cases can be collected in a plot shown in Fig. 1. In Fig. 1 we display in the xoj/o 2 o 
space the regions for the x 3 ± xy 2 germs and in Fig. 2 the separatrices among regions with different 
codimensions in parameter space. The separatrix xq = 1/2, say, in conjuction with the equations 
defining the instability surface (||) gives rise to the plane — ao + &o + ca = 2, used in Fig. 2 to find the 
separatrix in parameter space, which is a hyperbola branch. 




FIG. 1. Regions for x 3 — xy 2 potentials (inner triangle) and for x z + xy 2 potentials (3 outer triangles) 
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FIG. 2. Phase diagram showing the corank-2 instability surface. The intersection of the plane and the 
surface gives one of the three branches of the boundary between the x 3 + xy 2 and x 3 — xy 2 potentials. The 
diagonal corresponds to the 3-state Potts model. (See text) 

The equivalence between the 3-jet and the canonical germs implies that their codimensions are 
equal. Thus we have that 

c °d(j 3 (.9)) = cod{x 2 y) = oo, 
cod(j 3 (g)) = cod(x 3 -xy 2 )=3, 

cod(j 3 (<?)) = cod^ 3 + xy 2 ) = 3, (12) 

and we observe that the behavior of the codimension of the 3-jet is discontinuous according to Fig. 
1. Indeed, when the point (xo,yo) crosses the inner triangle the codimension of the 3-jet jumps to 
infinity. The codimension is only finite when what we could call hypothesis Si is fulfilled, that is, 

Ei = (l-2xo)(l-2|«,)(l-2^)^0. (13) 

On the contrary, if Si = 0, corresponding to the sides of the inner triangle, the germ g is not 3- 
determinate and we have to increase the order of its jet by one more degree to explore 4-determinacy. 
We have that 

fig) = aix 3 + a 2 x 2 y + a 3 xy 2 + a A y 3 + b x x A + b 2 x 3 y + b 3 x 2 y 2 + b A xy 3 + b 5 y 4 (14) 

where now the ^-coefficients are found to be: 

h = (x 3 + V 3 )/12, 

b 2 =z 3 /3, b 3 = z 3 /2, b 4 = z- 3 /3, 

b 5 = (y 3 + V 3 )/12. 

As before, in order to classify the 4-jet of g we invoke again a lemma of the Singularity Theory (see 
Appendix A). In all of the subcases examined the 4-jet is equivalent to the canonical germ x 2 y + y 4 , 
it is 4-dctcrminate and its codimension is equal to 4. 
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C. Determinacy and codimension of g(x,y) 

Once we have shown that the 3-jet is equivalent to the canonical germs x 3 ± xy 2 , and the 4-jet to 
x 2 y + y , it is clear that as these germs are 3- and 4-determinate, the 3-jet and 4-jet will be also 3- 
and 4-determinate, respectively. Then g ~ a; 3 — xy 2 (for Si < 0), g ~ x 3 + xy 2 (for Ei > 0) both 
with cod(.g) = 3, and g ~ x 2 y + y 4 (for Ei = 0) but with cod(.g) = 4. These equivalences exist up to 
an unknown change of coordinates, so that here and in the following the variables x and y have no 
physical interpretation and may be regarded as dummy variables. 

To see that x 3 ± xy 2 are 3-determinate and according to the main theorem about fc-determinacy, it 
suffices to prove that |l6|Jl7ll 

(x, y) 3+1 C (x, y) 2 A(x 3 ± xy 2 ) + (x, y) 3 + 2 (f 5) 

which is easy (for notation, see Appendix B). The same can be proven for the 4-determinate case but 
in Eq. ( [i"5| ) the exponent 3 should be replaced by 4. 



D. Canonical unfoldings of the germs x :i ± xy 2 and x 2 y + y 4 

It is well known from CT that {[x], [y], [x 2 ]} is a basis for the quotient vector spaces (x,y)/A(x 3 — 
xy 2 ) and (x,y)/A(x 3 +xy 2 ) and {[x], [y], [x 2 ], [y 2 ]} for (x,y)/A(x 2 y + y 4 ). Moreover, the /c-transversal 
(for all k > 0) canonical unfolding of the canonical form x 3 — xy 2 (elliptic umbilic, Z?_4) and its dual 
x 3 + xy 2 (hyperbolic umbilic, -D+4) and x 2 y + y 4 (parabolic umbilic, D 5 ) are respectively 

x 3 - xy 2 + Xix + \ 2 y + A 3 x 2 , 
x + xy 2 + Xix + \ 2 y + X3X 2 , 

and 

x2 y + y 4 + + ^y + A 3 x 2 + x^y 2 - 

The corresponding bifurcation diagrams are well known in the Singularity Theory fl8|| and can be 
seen in any of the standard books on this theory |^,^) . Equations governing such bifurcation diagrams 
are 

B D _ i : 3x 2 -y 2 + Xi+ 2X 3 x = 0, -2xy + X 2 = 0, 3x 2 + xX 3 + y 2 = 0, 
B D+i : 3x 2 + y 2 + Ai + 2X 3 x = 0, 2xy + A 2 = 0, 3.x 2 + xX 3 - y 2 = 0, 

and 

B D , : 2xy + Ai + 2X 3 x = 0, x 2 + 4j/ 3 + A 2 + 2X±y = 0, 6y 3 + yX A + QX 3 y 2 - x 2 + A3A4 = 0, 
which are obtained by equating to zero the first derivatives and Hessian of g in each case. 



E. Canonical 5-unfolding of g 

Since h has five parameters, the preceding unfoldings have to be extended with two or one irrelevant 
parameters in order to apply the isomorphy theorem. Thus, for example, for the canonical unfoldings 
x 3 ± xy 2 , the two new unfoldings denoted by /?i and (5 2 can be written as 
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0x(x, y, Ai, • • ■ , As) = x 3 - xy 2 + Xix + X 2 y + X 3 x 2 



(16) 



and 

fh(x, y, Ai, • • • , A 5 ) = x 3 + xy 2 + X^ + X 2 y + X 3 x 2 (17) 

and they are fc-transversal unfoldings for all k > of x 3 — xy 2 and x 3 + xy 2 , respectively. Now the 
new bifurcation sets are Bp 1 = Bp_ 4 x H 2 and Bp 2 = Br> +i x H 2 . 

In both cases, we can affirm that there is a change of coordinates (a diffeomorphism) <P( Xo ,u ) such 
that g = (x 3 — xy 2 ) ■ <P( Xo ,y ) an d the same holds for x 3 + xy 2 . Consequently, 

Pi(x,y,\ u ---,\ 5 ) = s 3 - st 2 + Ais + X 2 t + X 3 s 2 (18) 

and 

j3 2 (x,y,X 1 ,---,X 5 ) = s 3 + st 2 + X 1 s + X 2 t + X 3 s 2 , (19) 

where (s,t) = ip/ xo>yo -\(x,y), is a 3-tranversal unfolding of g with five parameters. Moreover, for 
the bifurcation sets we have that B^ = Bp 1 and B^ 2 = Bp 2 A similar reasoning can be used for the 
canonical unfolding x 2 y + y 4: + X\x + X 2 y + X 3 x 2 + A4J/ 2 but now we have only one irrelevant parameter. 

Finally, these canonical unfoldings of g need to be related to the function h or translated Gibbs 
potential. This is shown explicitely in next subsection through the so-called transversality condition. 



F. fc-transversality of the translated function h 



Now we return to the translated Gibbs potential h to establish its relation with the unfolding of 
g studied above. The condition for the existence of this relation is its transversality. Since we will 
be dealing with the physical function h its arguments will be the original physical variables, to be 
distiguished from the dummy variables used before. The vector space of the transversality, V/,, is 
defined by the first partial derivatives of h with respect to the five parameters according to 

V h = (D a h{x,y, 0,0, 0,0,0) - D a h(0,0, 0,0, 0,0,0), 
D b h(x, y, 0, 0, 0, 0, 0) - D b h(0, 0, 0, 0, 0, 0, 0), 
D c h(x, y, 0, 0, 0, 0, 0) - D c h(0, 0, 0, 0, 0, 0, 0), 
D^hix, y, 0, 0, 0, 0, 0) - D^hiO, 0, 0, 0, 0, 0, 0), 

D„ y h(x, y, 0, 0, 0, 0, 0) - D„ y h(0, 0, 0, 0, 0, 0, 0)> R (20) 

where (• • means all the linear combinations with real coefficients. 

Let us analyze the transversality of the unfolding h(x,y,a,b,c, n X i /-%)• As was mentioned in the 
Introduction, this study is carried out in order to show that the function h and the canonical unfolding 
of the germ g are isomorphic. The vector space of the transversality, Vh, is defined by Eq.(p0|) and 
reads in our case 

Vh = (-yx - y 2 + yz Q - y x - y Q y, -x 2 - xy + xz - x x - x y, xy + xy Q + x y, -x, -y)n- (21) 



Our next step is to prove that h is 3- or 4-transversal, according to each case. For this goal, we 
invoke one of the main theorems on fc-transversality which states that h will be fe-transversal when 
the following algebraic condition is met |jj6| , fi"7f 

(x,y) = A(g) + Vh + (x,y} k+1 , (22) 
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which is fulfilled in our case for k — 3 or 4. This requirement is proved in Appendix C with the 
hypothesis So and Si ^ or Si = 0, respectively. We finally conclude that ft, is a 3- or 4-transversal 
unfolding of g with five parameters. Moreover, by using the main theorem on fc-transversality, and 
for the 3-transversal case, (3\ and 02 & re isomorphic to h, that is, there are three diffeomorphisms and 
a perturbation of parameters depending on (xq, yo) for each case such that 

h(x, y, a, 6, c, /Lt x , n v ) = H(x + x ,y + yo, a + a , b + b a , c + c , fi x + ^,o, Ms/ + Mj/,o) 
- H (x , y ,a ,b ,c , fj, xfi , fi yfi ) 
= s 3 ± st 2 + Xis + X 2 t + A 3 s 2 + e( Xo ,y )(a, 6, c, MsiMy)> (23) 

where the diffeomorphisms take the following expressions: 

■^(x ,y )(x,y,a,b,c,n x ,fi y ) = (u,v,Xx,...,Xs), 
Vix^y^iu.v) = (s,i), 
V(x , Vo ){a,b,c,n x ,fi y ) = (Ai, . . . , A 5 ), (24) 

with 

iP (xo , yo) (x,y, 0,0, 0,0,0) = 0,0, 0,0,0). (25) 

All of these diffeomorphisms preserve the origin of coordinates. Moreover, concerning the bifurcations, 
we have that V(x ,yo)(Bh) = -Bd_ 4 x H 2 and r]^ Xo yo ^(Bh) = Bjj +l x M 2 . Similar expressions can be 
written for the 4-transversal case. 



IV. THE THREE-STATE POTTS MODEL 



The natural (and oldest) generalization of the Ising model consists of taking a site variable which can 
take three equivalent states instead of two, constituting the three-state Potts model. (Similarly, one 
can define the g-state Potts model.) It has complete permutation symmetry among the three states, 
yielding a Gibbs potential corresponding to (^) with a = b = c, Fig. 2, leaving as the only parameter the 
temperature T. The corank-2 instability occurs for T c = 1/9 = 0.1111 and the correponding potential 
is the D_4 germ. The 3-parameter unfolding ( |l6| ) contains the possible perturbations of temperature 
or chemical potentials. Since in the three-state Potts model only the temperature perturbation is 
allowed we must have a one-parameter unfolding. Its form is best deduced by symmetry arguments 
fiJI . We can substitute the germ by the symmetric form 

z 3 + z 3 = x 3 - ixy 2 7 (26) 

where we have introduced the complex variable z = x + i y. The permutation symmetry is obviously 
generated by the discrete rotations z — > e l s z and the reflection z — > z. The temperature perturbation 
preserves the symmetry and belongs to the vector space (x , y) / A(j 3 (g)) . Hence, it must be x 2 + y 2 — 
zz (see Appendix B). The corresponding unfolding is 

z 3 + z 3 + A z z = x 3 - 3 x y 2 + A (x 2 + y 2 ), (27) 

where A is a monotonic function of the temperature, at least, in a neighbourhood of T c . 

Now there arises the problem that the unfolding (pTj), being a section of the complete unfolding of 
-D_4, contains potentials with one minimum and three saddle points if A > or with one maximum and 
three saddle points if A < 0. These potentials are not bounded below nor have the three minima to be 
expected in this model. We must recall here the local character of CT and the discussion in |l| . There 
we remarked that the global character of these potentials can only be established by a numerical study 
over the entire range of the variables x and y. Thus one finds that they indeed have three minima , far 
from the point (0, 0) and distributed symmetrically, if we keep the temperature in a neighbourhood of 
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T c . For a large value of T only the minimum at (0, 0) survives. This is physically sensible, for at high 
temperature only the symmetric disordered phase can remain. The temperature T c precisely signals 
the point at which the symmetric disordered phase becomes unstable and disappears. Potentials for 
various situations are plotted in Fig. 3. 
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FIG. 3. Potentials for the 3-state Potts model at various temperatures. 



Note that for T > T c the three minima converge to (0,0), eventually merging there. Since their 
depth decreases as the depth of the minimum at (0, 0) increases, there must be a value of T such 
that they are equal. If we adopt the Maxwell convention then we have a first order phase transition. 
This phase transition is known to occur in this model in space dimension d — 4 and almost certainly 
for d > 3 p9[ . In d = 2 the phase transition is definitely of second order |12]| , which means that 
the mean-field theory potential with which we started is insufficient in low dimensions to describe 
accurately the real character of the possible phase transitions. 

Landau potentials for the 3-state Potts model that are bounded below have appeared in the literature 
p?j|,pl 12 1 . They are commonly derived by renormalization group (RG) arguments. In particular, it 
has been shown that for d > 3 the only RG relevant term further from the cubic is the quartic 
symmetric term (x 2 + y 2 ) 2 — (z z) 2 . It is even possible to prove with the powerful methods of 2d 
conformal field theory that this term, in addition to the cubic term, constitute a well defined Landau 



potential 12 1. So the potential c (z + z ) + [z z) seems to be adequate for this model. In fact, due to 



the addition of the quartic term, it exhibits the same behavior as the total non-polynomial potential 
displayed in Fig. 3. In any event, we must remark that the quartic term does not belong to the germ, 
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as calculated before. However, it could well arise from a more refined Gibbs potential, namely, of the 
type obtained with cluster variation methods B2 ], especially if symmetry arguments are invoked like 
in Ref. 0]. 



V. CONCLUSION 

We have analyzed the possible modes of instability and further singularities of the mean-field theory 
Gibbs potential for a three-component model in the case that two order parameters are needed for 
their description. Since the case of one order parameter has been already studied in Ref. JlJ], the 
study of this Gibbs potential is now complete. We have found that three new catastrophes take place, 
namely, the hyperbolic umbilic, Z?+4, its dual the elliptic umbilic, D-4, and the parabolic umbilic, D5, 
catastrophes. The parabolic umbilic has codimension 4, which is the highest we can reach in the case 
of two control parameters, unlike the case of one control parameter studied in jj], where codimension 
five was reached. The germ and unfolding of the elliptic umbilic catastrophe are precisely (isomorphic 
to) those of the 3-state Potts model, which is the model with the highest symmetry, and generally 
show three and even four phase coexistence. The total potential for the models belonging to the 
hyperbolic umbilic class do not have more than one local minimum and therefore they cannot give rise 
to phase coexistence. The phase structure for the parabolic umbilic is fairly complicated and can be 
seen in any of the standard books about CT |Hfl^8[. Finally, it is remarkable that with two order 
parameters the highest codimension 5, which should give rise to even more complex singularities, is 
not reached. 

We must say a few words about natural generalizations. Instead of limiting ourselves to three 
components we could well envisage the case of q components. This would include most models ever 
considered in the physics of phase transitions. The most symmetric case is the g-state Potts model. An 
essential advantage of increasing q is that mean-field theory becomes more accurate — in fact becoming 
exact for space dimension d > 2 if q > 4 [jHj. A Landau potential for the g-state Potts model with 
q — 1 order parameters is already known pi[ . Presumably, the rigorous study of the highest corank 
instabilities of the q component Gibbs potential would produce a germ equal to the cubic term of 
that Landau potential, as for the 3-state Potts model. The corresponding singularities exhibit the 
interesting feature of having modular parameters p3]| already for q — 4. There would certainly be a 
flock of other singularities with the same corank and the same or higher codimension, and with lower 
corank. We leave to the entrepreneurial reader the exploration of this endless world of mathematical 
entities and their physical application. 
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APPENDIX A: 



The application of the lemma of CT to Eq.(Jll|) to classify the 3-jet can be carried out in a systematic 
way considering the following subcases: 



(A) if 04 7^ (that is, 2y 7^ 1 — xo), then the auxiliary cubic equation of Eq.(|llJ), obtained by 
dehomogenizing the polynomial, is given by 

Q = Bl + ^ t + ^ + t s, (Al) 

O4 O4 O4 

Its discriminant is defined as A 2 = (ii — t 2 ) 2 (t 2 — tz) 2 (*3 — ^i) 2 > where t\, t 2 and t 3 are the roots of 
the cubic equation ([Al] ). It can be expressed as 

A 2 = -4p 3 r- 27 r 2 + 18pqr - 4q 3 + p 2 q 2 

= ^4 x o 4 % 4 *o~ 4 (1 - 2*o) (1 - 2yo) (1 - 2z ) (A2) 

with the following definitions: p = —03/04, q — 0,2/0,4 and r = —01/04. Notice the following subcases: 



(^4i) Eq.(Al) has three equal real roots (Aj = 0); this subcase is not possible because j/o > 0, 



(A 2 ) Eq.(Al) has three real roots but two of them are equal (that is, A 2 = 0); in this subcase, j 3 (g) 
and the monomial x 2 y are equivalent or, in mathematical terms, j 3 {g) ~ x 2 y (that is, the 3-jet and 
the monomial are equal up to a change of coordinates which essentially implies that 04 7^ 0), 

(A3) Eq. (|Al|) has three distinct real roots (that is, A 2 > 0) and j 3 (g) ~ x 3 — xy 2 , and 



(A4) Eq.(Al) has two conjugate complex roots (that is, Af < 0) and then j (g) ~ x 3 + xy , 



(B) if 04 = and a\ 7^ (that is, 2?/o = 1 — 2:0 and 2xq 7^ 1 — ya), then Eq.(Al) is replaced now by 
the auxiliary equation 

O^f+SV+t 3 (A3) 

Oi Oi 

and its discriminant A| is written as 

A 2 . = -4q 3 +p 2 q 2 

= K(x ,y )(2-3x -2y ) (A4) 

with K(xo,yo) > 0. The following redefinitions are now used: p = —02/0!, o = 03/02 ; £4, £ 2 and £3 
are again the new roots of the cubic equation, Eq. ([A3]). Thus we have again the following subcases 



(B\) Eq. (A3) has three equal real roots. This is not possible since 02 = 03 7^ 0, 



(B 2 ) Eq. (A3) has three real roots, two of them equal (A 2 = 0); this condition implies that 
3x + 2y = 2 and then j 3 (g) ~ x 2 y, 

(B 3 ) Eq. (|A3|) has three distinct real roots (A| > 0), then 2 > 3x + 2y and j 3 (g) ~ x 3 - xy 2 , and 



(B4) Eq. (A3) has two conjugate complex roots (A 2 , < 0), then 2 < 3a;o + 2yo and j 3 (<?) ~ x 3 +xy 2 
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(C) Finally we have the case where 04 = a\ = (or xq — yo = 1/3) and we obtain that j 3 (g) ~ x 3 —xy 2 . 
Due to the fact this case is related to the well known Potts model, a more detailed analysis of this 
case will be addressed in Section IV. 



For the 4-jet (Si = 0), the following subcases can be considered: 

(Di) if 04 7^ (that is, 2yo ^ 1 — xq), then we take the auxiliary cubic equation 

= at + a 2 t + a 3 t 2 + a 4 t 3 . (A5) 

Eq. ( |AE| ) can have only three real roots of which two are equal. With the linear transformation 

ip = u = —a,4 t\ x + 04 y, v — — t 2 X + y (A6) 

with t\ the single root and t 2 the double root of Eq. ( |A5| ) , then 

f(g)i;(u,v) =uv 2 . (A7) 

The value for the double root t 2 is easily obtained from Eq. ( |A5| ) and its derivative with respect to 
the variable t since both equations are satisfied for t 2 \ 

, „ . . „ ( -2 w if x = 1/2 

a 2 a 3 -9aiffl 4 1 + 2 .t (-1 + y ) - 2 y I l if „ _ i /o fA o^ 

2(3a 2 a 4 -a 3 2 ) 2x 2 [ _i + _i_ if Zq = i /2 

One can further obtain t\ and it coincides with t 2 only on the vertices of the triangle, where a 
concentration vanishes. Therefore, this situation, which would lead to a different singularity with 
3-jet u 3 , is actually outside the domain we consider. 

We can write the 4-jet in the new variables as 

j 4 (g) i>(u, v) = uv 2 + ciu 4 + C2U 3 v + czu 2 v 2 + auv 3 + c^v 4 , (A9) 

where the new coefficients Ci are linear combinations of the bi. Only the sign of the non-zero c\ 
coefficient is needed since, according to the mentioned lemma, the 4-jet is equivalent to 

f(g) ~x 2 y + sign( Cl )y 4 , (A10) 



provided that c± ^ 0. This comes from the fact that all the 4-degree monomials of the 4-jet ( A9) belong 
to the Jacobian ideal of u v 2 except precisely u and therefore can be removed by a diffeomorphism. 
From the linear transformation Eq. ( |A6| ), we have that 

ci = 7 4 [h + b 2 t 2 + b 3 t 2 + bit 3 + b 5 t 4 ] (All) 



with 7 = l/[a,4(t 2 —ti)]. Therefore the sign of c\ is given by the factor inside the bracket. Substituting 



for t 2 , c\ = 7 4 x 4 and is always positive, for any of the three factors given by Si = 0, that is: yo = 1/2, 



Xq = 1/2 or zq = 1/2. Finally, 

cod(j 4 ( 5 )) = cod(x 2 y + y 4 ) = 4 (A12) 

(D 2 ) if a4 = 0, ai ^ (that is, 2yo = 1 — xq, 2xq ^ 1 — yo) and Si = 0, then only one point (xq, yo) 
needs to be studied, Xo = 1/2 and yo = 1/4. The corresponding auxiliary cubic equation is now 

= t 3 + ^t 2 + ^t. (A13) 
ai ai 
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The roots of this last equation are: t\ = and t<i — —2. With the linear transformation 

vjr 1 = u = x + 2y, v = 2x (A14) 

applied to the 4-jet we obtain that 

j 4 (g) ip(u, v) = vu 2 + rfiw 4 + d 2 u 3 v + d 3 u 2 v 2 + d 4 uv 3 + d 5 v A . (A15) 

Here the sign of d§ is positive and therefore the 4-jet is again equivalent to the canonical germ x 2 y+y , 
its codimension being again equal to 4. 



APPENDIX B: 

In this Appendix we are going to present the mathematical concepts and notation widely introduced 
m Refs. yjlj and necessary to follow the main steps developed in Section III. 

Let us consider real functions of class oo and defined in a neighbourhood of £ H™. We establish 
that two functions are equivalent if they coincide in a neighbourhood of 0. The classes we obtain 
are called germs of functions and the set of germs is denoted by E(n). The operations f + g and 
f g give to E(n) the estructure of a ring and M(n) = {/ G E(n)/f(0) — 0} is a maximal ideal 
of this ring. Moreover, the operations f + g and A/ with A £ IR give to E(n) the structure of 
a real vector space of dimension oo. The ideal M(n) is generated by x\, . . . ,x n , that is, M(n) = 

{fixi H h f n x n / /i, . . . , /„ £ E(n)}. In general, if /i, . . . , /„ € E(n), we designate by . . . , /„} 

the ideal generated by f±, . . . , /„, that is, 

■■■,/«) = {/iSfiH \-fn9n /gu---,9n eE(n)}. 

In particular, M(n) = (x\, . . . ,x n ). 

It is possible to define powers of M{n) as M(n) k . It can be proven that M(n) k is equal to the 
ideal of E(n) generated by the monomials in Xi,...,x n of degree k. In particular, for example, 
(x,y) 2 = (x 2 ,xy,y 2 ) and 

(x, y) 2 = {f x x 2 + f 2 xy + f 3 y 2 / f 1} f 2 , h G E{2)}. (Bl) 

We have also that M(n) fe+1 = {/ € E(n) / D l f(0) = 0, i < k} where by D l we mean the derivative 
of degree i. 

Now the main theorem about fc-dctcrminacy establishes a sufficient condition for a germ g to be 
/c-determinate. This condition reads 

(x,y) k+1 C (x,y) 2 .A(g) + (x,y) k + 2 , (B2) 

and, in our case, it will be applied for k — 3 and 4. The factor A (g) is known as the ideal of Jacobi and 
is defined from the first partial derivatives of g with respect to x and y (D x g and D y g, respectively). 
The calculation of this ideal is the general starting point to establish the determinacy and codimension 
of g. Thus, for example, 

A (g) = (D x g(x,y),D y g(x,y)) 

= (^x 2 M + xyz Q 2 + \y 2 zv 2 + \^{x^ 3 + zo 3 ) + x 2 yz^ 3 + xy 2 z^ 3 + |y 3 ^ 3 + r(x, y), 

i x 2 z 2 + xy z ~ 2 + ^y 2 N + \x 3 z^ 3 + x 2 yz^ 3 + xy 2 z^ 3 + ^y 3 (% 3 + z ~ 3 ) + V)), ( B3 ) 

with the definitions M = —x^ 2 + z^ 2 and = —y^ 2 + Zq 2 , and with r(x,y) , s(x,y) £ (x, j/) 4 after 
Taylor expansions of D x g and D y g around the point (0, 0) have been performed, the first terms being 
ignored according to the hypothesis So- Analogously, we have that 
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A(x 3 ±xy 2 ) = (3x 2 ±y 2 ,±2xy). (B4) 

One is not always in the position to invoke a given lemma from CT to classify the germ and establish 
its fc-d eter minacy In what follows we illustrate how to proceed in the general case using the condition 



above (B2). Thus, 



(x, y) 2 A(g) + (x, y) 5 = (x, y) 2 {^x 2 M + xyz^ 2 + ~V 2 Zq 2 , 
\ x 2 Zq 2 + xyz 2 + ^y 2 N) + (x,y) 5 
= {x,V? + {\x*M + x 3 yz^ 2 + \ l x 2 y 2 z^ 2 1 

- x 4 z~ 2 + x 3 yz- 2 + ^x 2 y 2 N, 

- x 2 y 2 M + xy 3 z- 2 + \v A Zq 2 , ^x 2 y 2 z~ 2 + xy 3 z~ 2 + ^y 4 N, 

- x 3 yM + x 2 y 2 z 2 + ^xy 3 z~ 2 , 

- x 3 yz~ 2 + x 2 y 2 z- 2 + ^xy 3 N), (B5) 
and therefore we have to prove that 

Hix 4 + fi2X 3 y + [i?,x 2 y 2 + fi4xy 3 + ^y* = vi[-x A M + x 3 y z^ 2 + -x 2 y 2 z^ 2 ] 

+ M\ x4>z o 2 + x3 V z o 2 + 7^x 2 y 2 N] 

i r ^ 2 2 71 , 3-2,1-4 -2i 

+ m-^x y M + xy a z + -y z } 
+ v^x 2 y 2 z- 2 + xy 3 z 2 + ^N] 



+■ v 5 {^x 3 yM + x 2 y 2 z 2 + ^xy 3 z 2 ] 



+ v 6 [lx 3 y z 2 + x 2 y 2 z 2 + ^xy 3 N] (B6) 

and by equating coefficients of the same degree we obtain five equations relating the \i coefficients with 
the v coefficients. This linear system of equations has solutions in Ui, . . . , v§ only if the hypothesis Si 
is fulfilled, Eq.@. 

Under these conditions, the germ and the 3-jet of the germ are equivalent, g ~ j (g) and cod(<?) = 
cod(j 3 (g)) = dim. vect. (x , y) / A(j 3 (g)) , that is, the dimension of the quotient vector space associated 
by the ideal of Jacobi of the 3-jet. Moreover, it can be also proved that {[x], [y], [x 2 ]} is a basis of that 
quotient vector space. Similarly, alternative basis can also be {[x], [y], [y 2 ]} or {[x], [y], [x 2 + y 2 ]}. 

The same procedure can be applied to the 4-determinate case. 



APPENDIX C: 



As has been mentioned above we need to show that h is 3-transversal (similar calculations are 
needed in order to show that for the hypothesis Si = 0, h is 4-transversal) . One of the main theorems 
about transversality establishes that this property is fulfilled when 

(x,y)=A(g) + V h + {x,y) 3+1 . (CI) 
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The ideal of Jacobi of g is given by Eq. (|B3j) and the vector space Vh by Eq.(|21]). 

In Eq.(|B3]), the ideal of Jacobi of g contains monomials of degree greater than 2. We are going to 
show that monomials of degree equal to 3 belong to A(j 3 (g)). Thus we have that 

Aix 3 + X 2 x 2 y + X 3 xy 2 + X 4 y 3 

= (ax + by)(^x 2 M + xyz^ 2 + \y 2 ^ 2 ) + [ex + dy)(^x 2 z^ 2 + xyz^ 2 + \y 2 N) (C2) 

and the linear system of equations obtained from equating coefficients of the same degree has solutions 
in a, b, c and d if Ei ^ . 

On the other hand, a straightforward consequence of what we have shown above is that 

Mg) + (x,y) 4 = A(f(g)) + (x,y)\ (C3) 
So finally we can rewrite the transversality condition as 

( x > V) = (^M + xyz^ 2 + 7^y 2 Zo 2 , \ x 2 Zq 2 + xyz Q 2 + ^ y 2 N) + (x, y} 4 

+ {~yx - y 2 + yz - yox - y y, ~x 2 - xy + xz - x x - x y, xy + xy + x y, -x, -y)n- (C4) 

In other words, we have to find a set of ^-parameters fulfilling 

[HX + [i 2 y + Hzx 2 + /i 4 y 2 + [i 5 xy = v x (-yx - y 2 + yz Q - y x - y y) 

+ v 2 {-x 2 - xy + xz Q - x x - x y) 
+ v z (xy + xy + x y) 
+ Vi(—X) 

+ vs{-y) 

+ "6 [^x 2 M + xyo^ 2 + -y 2 a^ 2 } 

+ v 7 [^x 2 ao 2 + xya^ 2 + \y 2 N]. (C5) 

Again by equating coefficients of the same degree we obtain a set of equations between the family of 
known /i and unknown v parameters. The corresponding system of equations has a matrix of rank 
equal to 5 and the transversality condition is fulfilled. 

For the 4-transversal case, the procedure is entirely similar. Only the following observation needs 
to be made. Monomials of degree 1 and 2 are included in the vector space of transversality and 
monomials of degree 5 and higher are considered in the term (x,y) 5 . So the equivalent to Eq. (C5) 
has to take into account monomials of degree 3 and 4. This leads to a system of nine equations and 
ten unknown parameters. The rank of that system is 9. 
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